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Abstract.  We  present  the  first  maximum-entropy  maps  re¬ 
constructed  from  visibility  amplitudes  obtained  with  an  opti¬ 
cal  long-baseline  interferometer.  They  show  the  Be  shell  star 
£  Tauri.  In  the  continuum  at  550  nm  the  star  is  an  unresolved 
point  source,  whereas  in  the  light  of  the  Ho  line  an  extended, 
elongated  structure  is  observed.  The  best-fit  parameters  of  an 
elliptical  Gaussian  model  are:  major  axis  a  =  3.55  ±0.33  mas, 
axial  ratio  r  =  0.30  ±  0.03,  position  angle  <j>  =  —59°  ±  4°. 
The  morphology  of  the  Ha  emission  region  is  most  easily  in¬ 
terpreted  as  a  disk  seen  almost  edge-on;  if  this  disk  is  axially 
symmetric,  a  limit  for  the  inclination  i  >  73°  can  be  derived. 
Our  observations  provide  further  evidence  for  the  equatorial 
disk  model  of  Be  stars. 

Key  words:  Techniques:  interferometric  -  Stars:  emission  line, 
Be  -  Stars:  imaging  -  Stars:  individual:  C  Tau  -  Stars:  mass-loss 


1.  Introduction 

While  stellar  long-baseline  interferometry  has  become  a  very 
successful  tool  to  measure  stellar  diameters  (e.g.  Mozurkewich 
et  al.  1991)  and  to  determine  “visual”  orbits  of  close  binaries 
(e.g.  Armstrong  1992),  extraction  of  more  detailed  informa¬ 
tion  about  the  stellar  brightness  distribution  from  the  visibil¬ 
ity  data  remains  a  considerable  challenge.  The  first  notable 
achievement  in  this  area  was  the  measurement  of  Sirius’  limb 
darkening  by  Hanbury  Brown  et  al.  (1974)  with  the  Narrabri 
intensity  interferometer.  Possible  structure  in  f)  Andromedae 
was  inferred  from  I2T  data  at  1.65  and  2.2 /rm  by  Di  Benedetto 
and  Bonneau  (1990).  Observations  with  the  Mklll  interferome¬ 
ter  have  shown  that  Mira’s  atmosphere  is  strongly  asymmetric 
and  changes  shape  on  the  timescale  of  months  (Quirrenbach  et 
al.  1992).  The  Ha  emission  region  of  the  Be  star  7  Cassiopeiae 
has  been  observed  with  the  I2T  (Thom  et  al.  1986),  GI2T 
(Mourard  et  al.  1989),  and  Mklll  (Quirrenbach  et  al.  1993) 
instruments;  the  results  support  the  long-standing  picture  of  a 
rotating  disk-shaped  envelope. 
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Whereas  the  interpretation  of  these  data  sets  has  been 
based  on  fitting  models  with  a  small  number  of  free  param¬ 
eters  to  the  measured  visibilities,  future  interferometric  arrays 
will  rely  on  methods  similar  to  those  currently  used  in  radio 
interferometry  to  reconstruct  the  sky  brightness  distribution. 
These  algorithms  have  already  been  adapted  for  the  analysis  of 
data  obtained  with  aperture  masks  on  large  telescopes  (Haniff 
et  al.  1988,  Buscher  et  al.  1990).  Closure  phase  information 
will  be  crucial  to  obtain  images  with  high  fidelity  and  dynamic 
range  using  instruments  such  as  the  Big  Optical  Array1  (BOA), 
the  ESO  Very  Large  Telescope  Interferometer  (VLTI)  and  the 
Cambridge  Optical  Aperture  Synthesis  Telescope  (COAST).  It 
is  possible,  however,  to  obtain  images  of  simple  objects  from 
amplitude  data  only.  Here  we  present  a  maximum-entropy  re¬ 
construction  of  the  Ha  disk  of  the  Be  shell  star  (  Tauri  from 
data  obtained  with  the  Mklll  Optical  Interferometer. 

2.  The  observations 

C  Tauri  (B4  Hie,  my  =  3.0)  was  observed  with  the  Mklll 
instrument2  on  12  nights  between  Sept.  27  and  Nov.  30,  1992; 
some  test  observations  had  been  performed  on  Dec.  27,  1991, 
and  Jan.  01,  1992.  The  Mklll  is  a  Michelson  interferometer 
with  a  single  north-south  baseline,  whose  length  can  be  varied 
between  3.0  and  31.5  m.  In  addition  to  the  wide-band  channel 
used  for  active  fringe-tracking,  three  data  channels  are  avail¬ 
able;  their  bandpasses  are  defined  by  interference  filters.  Auto¬ 
matic  acquisition  of  the  stellar  images  and  of  the  fringes  under 
computer  control  allows  for  rapid  switching  between  different 
stars,  so  that  observations  of  the  program  stars  can  be  sand¬ 
wiched  between  calibrators  with  known  diameters.  With  this 
observing  strategy,  an  accurate  calibration  of  the  system  visi¬ 
bility  (i.e.,  the  fringe  contrast  obtained  on  a  point  source),  can 
be  achieved  as  a  function  of  seeing,  zenith  angle,  and  time.  De¬ 
tailed  descriptions  of  the  instrument  and  observing  procedures 
are  given  by  Shao  et  al.  (1988)  and  Mozurkewich  et  al.  (1991). 

'The  BOA  is  the  imaging  portion  of  the  Navy  Prototype  Op¬ 
tical  Interferometer  under  construction  at  Lowell  Observatory. 

2The  Mklll  Optical  Interferometer  on  Mt.  Wilson  near  Los 
Angeles,  CA,  is  operated  by  the  Remote  Sensing  Division  of 
the  Naval  Research  Laboratory  (NRL). 
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For  the  observations  reported  here,  we  used  baselines  rang- 
^  ing  from  4.2  to  31.5  m.  The  main  data  channel  was  centered 
^  on  the  Ha  line  at  656  nm  and  had  a  width  of  1  nm.  A  second 
cri  channel  at  the  same  central  wavelength,  but  10  nm  wide,  was 
:  1  used  for  control  purposes;  this  is  necessary  because  the  1  nm. 
wide  filter  tends  to  attenuate  the  wings  of  the  rather  wide 
Ha  line.  The  third  channel  was  centered  at  550  nm  and  25  nm 
wide;  it  was  used  to  measure  the  stellar  continuum  for  com¬ 
parison.  Each  individual  observation  (“scan”)  was  300s  long; 
this  gave  sufficient  signal-to-noise  on  the  program  as  well  as 
calibrator  stars  in  all  channels.  The  observations  were  sched¬ 
uled  in  a  way  that  would  optimize  the  hour-angle  coverage  on 
the  longer  baselines;  the  resulting  coverage  of  the  Fourier  (uv) 
plane  for  the  combined  data  on  £  Tauri  from  all  nights  in  the 
narrow  Ha  channel  is  shown  in  Fig.  1. 
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Fig.  1.  Coverage  of  the  Fourier  (uv)  plane,  for  the  data  in  the  Ha 
channel. 


3.  Maximum-entropy  mapping 

The  Mklll  standard  data  reduction  software  incoherently  in¬ 
tegrates  the  square  of  the  visibility  amplitude  V^v,  estimated 
from  each  4  ms  sample  over  the  length  of  the  scan,  and  divides 
it  by  the  system  visibility,  giving  V^\ .  The  calibrated  data  can 
then  be  used  by  model  fitting  or  image  reconstruction  programs 
similar  to  those  used  in  radio  interferometry.  For  the  analysis 
of  the  £  Tauri  data,  we  used  an  iterative  maximum-entropy 
algorithm,  which  seeks  to  find  the  “smoothest”  image  which 
is  compatible  with  the  data.  In  each  step,  it  calculates  model 
visibilities  from  a  test  image,  compares  them  to  the  data  and 
computes  and  its  partial  derivatives,  and  modifies  the  im¬ 
age  along  the  direction  of  the  gradient  of  Xied  ~  AS,  where  Xicd 
is  the  reduced  x2i  A  a  Lagrange  multiplier,  and  S  the  image 
entropy  (see  e.g.  Narayan  and  Nitnayanda  1984).  Starting  from 
a  user-specified  initial  image  (typically  a  point  source  at  the 
origin  on  a  low  uniform  background),  the  algorithm  thus  tries 
to  converge  to  the  image  which  has  y(ed  =  1  and  maximum 
entropy.  More  detailed  descriptions  of  the  algorithm  are  given 
by  Skilling  and  Bryan  (1984)  and  Buscher  (1993). 
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Fig.  2.  Maximum-entropy  reconstruction  of  £  Tauri  in  the  contin¬ 
uum  channel  centered  at  550 nm.  Contour  levels  are  0.1,  0.2,  0.5,  1, 
2,  5,  10,  20,  50,  and  80%  of  the  peak. 


milliarcseconds 

Fig.  3.  Maximum-entropy  reconstruction  of  £  Tauri  in  the  Ha  emis¬ 
sion  line  (656nm).  Contour  levels  are  as  in  Fig.  2.  The  arrow  indi¬ 
cates  the  position  angle  of  the  linear  polarization. 


We  tested  the  program  using  one  night’s  data  on  the  bi¬ 
nary  star  Capella  (a  Aurigae);  it  converged  to  an  image  with 
two  pointlike  stars  at  a  separation  and  position  angle  in  agree¬ 
ment  with  results  from  model  fitting  and  from  a  new  improved 
orbital  solution  (Hummel  et  al.  1994).  During  these  tests  we 
found  that  the  convergence  was  substantially  improved  and  the 
sidelobe  level  in  the  final  image  drastically  reduced  when  V 2 
at  zero  spacing  was  explicitly  set  to  1.  In  principle,  this  effect 
is  well-known  in  radio  astronomical  applications,  but  there  the 
flux  density  at  zero  spacing  is  frequently  not  known  precisely 
enough  to  aid  the  maximum-entropy  algorithm  significantly. 

The  photosphere  of  £  Tauri  is  expected  to  be  almost  un- 
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“  resolved  even  on  our  longest  baseline;  estimates  for  the  stel¬ 
lar  angular  diameter,  based  on  photometric  data,  range  from 
0.4  mas  (Ochsenbein  and  Halbwachs  1982)  to  0.6  mas  (Glush- 
«j  neva  1987).  Our  data  in  the  550  nm  channel  show  indeed 
g)  only  marginal  indications  of  resolution;  on  the  31.5  m  base- 
'  1  line  V2  «  0.8.  This  allows  us  to  set  an  upper  limit  of  r;  1  mas 
on  the  size  of  the  photosphere.  Based  on  radial  velocity  varia¬ 
tions  of  the  shell  lines  it  has  been  argued  that  £  Tauri  is  in  a 
binary  system  with  period  132.91  days  (Delplace  and  Cham- 
bon  1976).  Since  there  is  no  signature  of  the  secondary  in  our 
550  nm  data,  we  can  set  a  limit  on  the  magnitude  difference 
A mv  >  2.5  (unless  the  separation  d  >  100  mas,  which  is  in¬ 
compatible  with  the  total  mass,  period  and  distance  of  the 
system).  This  is  consistent  with  the  finding  by  Floquet  et  al. 
(1989)  that  the  presumed  secondary  cannot  be  a  bright  (lumi¬ 
nosity  class  II)  giant,  which  would  fill  its  Roche  lobe. 

In  Fig.  2,  we  present  the  maximum-entropy  map  of  £  Tauri 
at  550  nm.  This  map  is  basically  the  point-source  response  cor¬ 
responding  to  our  uv  coverage;  the  apparent  elongation  in  E-W 
direction  is  due  to  the  fact  that  the  resolution  is  substantially 
better  in  the  N-S  than  in  E-W  direction  (see  Fig.  1). 

In  the  light  of  the  Ha  line,  £  Tauri  is  clearly  resolved  and 
strongly  elongated  (in  the  projection  onto  the  sky).  Figure  3 
shows  the  maximum-entropy  map  from  the  data  taken  in  the 
1  nm  wide  channel.  The  fit  to  the  data  is  excellent,  as  can 
be  seen  in  Fig.  4,  where  we  superpose  the  model  prediction 
on  the  data  obtained  on  the  longest  (31.5  m)  baseline.  (The 
baseline  position  angle  is  the  angle  of  the  corresponding  point 
in  the  uv  plane,  measured  clockwise  from  the  v  axis.)  It  should 
be  pointed  out  that  the  strong  elongation  of  the  Ha  emission 
region  is  apparent  directly  from  the  visibility  data:  a  spherically 
symmetric  object  would  give  visibilities  that  are  symmetric 
with  respect  to  meridian  transit  (baseline  position  angle  0°). 


Fig.  4.  Visibility  data  in  the  1  nm  wide  Ha  channel  from  the  longest 
(31.5m)  baseline,  and  maximum-entropy  model. 

We  have  performed  a  series  of  tests  with  a  number  of  dif¬ 
ferent  starting  models  to  verify  that  the  maximum-entropy  al¬ 
gorithm  converged  consistently  to  basically  identical  solutions. 
Nevertheless,  it  should  be  borne  in  mind  that  a  maximum- 


entropy  “image”  is  nothing  but  a  model  representation  for  the 
data  and  thus  susceptible  to  a  number  of  artifacts. 

First,  the  lack  of  phase  information  in  our  data  will  tend 
to  favor  a  mirror-symmetric  reconstruction.  Our  test  with  a 
reconstruction  of  the  binary  star  Capella,  which  did  not  give  a 
“ghost  image”  of  the  secondary,  shows  that  this  is  not  necessar¬ 
ily  the  case,  however.  If  the  object  is  fairly  simple,  amplitude 
data  with  sufficient  signal-to-noise  can  produce  asymmetric  re¬ 
constructions.  In  the  case  of  £  Tauri  we  cannot  decide  whether 
the  structure  of  the  envelope  is  indeed  symmetric  or  whether 
the  signal-to-noise  in  the  data  is  not  high  enough  to  make  the 
maximum-entropy  program  converge  to  an  asymmetric  solu¬ 
tion.  Rapid  variations  of  the  emission  line  profiles  (e.g.  Yang 
et  al.  1990)  suggest  that  the  Ha  emitting  region  might  be  ei¬ 
ther  asymmetric  or  not  centered  on  the  star  (Telting  et  al. 
1993).  On  the  other  hand,  Thom  et  al.  (1986)  find  that  the 
envelope  of  7  Cassiopeiae  is  in  fact  symmetric  and  centered  on 
the  star  to  within  0.3  mas.  Future  interferometric  observations 
with  longer  baselines  and  high-  quality  phase  information  are 
needed  to  decide  this  question. 

Second,  the  asymmetric  coverage  of  the  uv  plane  will  show 
up  as  an  artificial  elongation  of  the  contour  lines  in  E-W  direc¬ 
tion,  as  discussed  above.  The  continuum  map  in  Fig.  2  can  be 
used  as  a  reference  for  the  point-source  response.  It  should  be 
pointed  out  that  maximum-entropy  maps  do  not  have  a  defined 
“beam  shape”.  (To  facilitate  the  comparison,  the  reproductions 
of  the  two  maps  have  been  scaled  by  the  wavelength  ratio.) 

To  derive  numerical  parameters  for  the  Ha  emission  region, 
we  have  therefore  also  fitted  a  simple  elliptical  Gaussian  model 
to  the  visibility  data.  The  best-fit  parameters  of  this  model 
are:  major  axis  a  =  3.55  ±  0.33  mas,  axial  ratio  r  =  0.30  ±0.03, 
position  angle  1 j>  —  —59°  ±  4°  (measured  from  north  through 
east).  The  quoted  errors  are  formal  errors  obtained  from  the 
increase  of  x2  when  the  parameters  are  varied;  the  value  for  r 
might  actually  be  smaller,  since  the  size  of  the  minor  axis  is 
very  close  to  our  resolution  limit.  This  elliptical  model  gives 
a  very  good  representation  of  the  overall  appearance  of  the 
maximum  entropy  map. 

The  total  extent  of  the  Ha  emission  region  is  «  10  mas  (3 
FWHM  of  the  Gaussian  model,  see  also  Fig.  3),  or  about  15 
photospheric  diameters.  This  value  agrees  roughly  with  model 
predictions  (e.g.  Poeckert  and  Marlborough  1978). 

4.  Discussion 

The  observed  axial  ratio  of  the  Ha  emission  region  r  =  0.30 
implies  that  the  intrinsic  geometry  is  strongly  non-spherical. 
This  result  is  inconsistent  with  nearly  spherical  models  for  Be 
stars  (e.g.  Doazan  1987)  and  provides  additional  evidence  for 
a  disklike  structure.  Both  fairly  thick  (opening  angle  ~  15°, 
Poeckert  and  Marlborough  1978)  and  very  thin  (opening  angle 
~  0.5°,  Bjorkman  and  Cassinelli  1993)  disks  are  consistent  with 
our  data;  in  the  latter  model  it  might  be  difficult  to  explain 
the  strong  shell  absorption,  however. 

Under  the  assumption  that  the  disk  is  axially  symmetric,  a 
lower  limit  to  the  inclination  of  the  axis  i  >  arccos0.30  =  73° 
is  derived  from  a  simple  geometric  consideration.  The  high  in¬ 
trinsic  degree  of  polarization  (Poeckert  et  al.  1979,  McDavid 
1990)  provides  independent,  although  less  direct,  evidence  for 
a  large  inclination  of  £  Tauri.  Published  values  for  v  sin  i  range 
from  220kms_1  (Slettebak  1982)  to  320kms_1  (Uesugi  and 
Fukuda  1982),  but  the  upper  end  of  this  range  seems  to  be 
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preferred  (Yang  et  al.  1990).  Gao  and  Cao  (1986)  have  tried  to 
determine  the  inclination  of  £  Tauri  using  a  method  proposed 
by  Hutchings  (1976),  which  is  based  on  the  comparison  of  op¬ 
tical  and  UV  line  widths;  they  found  i  60°.  Combining  these 
values  for  v  sin  i  and  i  leads  to  conclusion  that  the  rotation 
rate  of  (  Tauri  is  very  close  to  the  critical  value,  a>/u>crit  =  1 
(Ruusalepp  1989).  The  larger  inclination  indicated  by  our  data 
gives  w/wcrit  «  0.85,  however,  supporting  the  idea  that  Be  stars 
do  not  rotate  at  the  critical  velocity.  Our  limit  for  the  inclina¬ 
tion  is  also  consistent  with  the  supposition  that  Be  stars  with 
shell  absorption  lines  in  the  visible  part  of  the  spectrum  are 
seen  almost  equator-on  (e.g.  Slettebak  1979). 

The  continuum  polarization  of  Be  stars  is  believed  to  be 
due  to  scattering  by  free  electrons  in  the  circumstellar  disk 
(McLean  1979).  This  explanation  predicts  that  the  polariza¬ 
tion  position  angle  x  should  be  perpendicular  to  the  equatorial 
plane.  Poeckert  et  al.  (1979)  and  McDavid  (1990)  find  x  ~  33°; 
with  only  a  small  interstellar  contribution.  The  position  angle 
<t>  =  ~59°  ±  4°  from  our  data  is  indeed  perpendicular  to  \- 

5.  Conclusions 

For  the  first  time,  we  have  been  able  to  reconstruct  a  max¬ 
imum-entropy  map  from  visibility  amplitudes  obtained  with 
optical  long-baseline  interferometry.  This  represents  an  impor¬ 
tant  step  towards  routine  imaging  with  milliarcsecond  resolu¬ 
tion,  the  declared  goal  of  the  interferometric  arrays  currently 
under  construction. 

Our  maximum-entropy  map  of  (  Tauri  in  the  Ho  line 
should  be  interpreted  with  some  caution,  because  of  the  lack 
of  phase  information  and  the  non-uniform  coverage  of  the  uv 
plane.  Nevertheless,  it  confirms  several  predictions  of  the  equa¬ 
torial  disk  model  for  Be  stars.  While  elongated  structures  have 
been  observed  for  the  radio  emission  of  ip  Persei  (Dougherty 
and  Taylor  1992)  and  for  the  Ha  emission  of  7  Cassiopeiae 
(Quirrenbach  et  al.  1993),  the  observed  axial  ratio  is  more  ex¬ 
treme  for  (  Tauri.  This  allows  tighter  constraints  to  be  placed 
on  the  intrinsic  geometry;  a  meaningful  lower  limit  to  the  incli¬ 
nation  can  be  derived  under  the  assumption  of  axial  symmetry. 

Future  observations  with  longer  baselines  and  improved  uv 
coverage  will  give  valuable  additional  information  about'  the 
shape  and  structure  of  Be  star  disks.  The  emission  line  regions 
of  these  and  similar  objects  are  particularly  good  targets  of  op¬ 
tical  long-baseline  arrays.  The  small  central  star  can  be  used  as 
a  reference  object  for  fringe  tracking,  and  the  line— continuum 
phase  difference  can  be  used  directly  for  imaging;  phase  closure 
methods  are  not  needed.  This  gives  a  large  increase  in  phase  in¬ 
formation  for  arrays  with  a  small  number  of  elements.  Further¬ 
more,  even  data  with  moderate  signal-to-noise  will  carry  valu¬ 
able  astrophysical  information.  (In  contrast,  the  main  goal  of 
the  BOA,  namely  imaging  of  stellar  photospheres,  requires  ex¬ 
tremely  high  signal-to-noise  and  excellent  calibration  to  detect 
small  deviations  from  simple  uniform  or  limb-darkened  disk 
models.)  Finally,  a  simple  estimate  shows  that  the  BOA  will 
be  sensitive  enough  to  put  ~  10  spectral  channels  across  the 
Ha  line.  A  combination  of  the  approaches  taken  by  Mourard 
et  al.  (1989)  and  in  this  work  will  thus  make  detailed  studies 
of  the  kinematics  of  Be  star  disks  possible. 
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